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Abstract 

We review the current high-significance X-ray detections of Warm- Hot Intergalactic 
Medium (WHIM) filaments at ^; > along the lines of sight to the two blazars 
Mrk 421 {z = 0.03) and lES 1028+511 {z = 0.361). For these WHIM filaments, we 
derive ionization corrections and, when possible, metallicity estimates. This allows 
us to obtain refined estimates of the number density of O VII WHIM systems down 
to the O VII column density sensitivity of our observations, and most importantly, 
a measurement of the cosmological mass density ri^^^'^ in the WHIM, at redshift 
z < 0.361. These estimates agree well with model predictions and with the total 
estimated amount of missing baryons in the local Universe, although errors are 
large, due to the still limited number of systems. We conclude discussing future 
observational strategies and mission designs for WHIM studies. 
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1 'Missing Baryons' and the Warm-Hot Intergalactic Medium 



The total number of baryons in the Universe, as inferred by both Big-Bang Nu- 
cleosynthesis (BBN) compared with observations of light element ratios (Kirk- 
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man ct al., 2003) and Cosmic Microwave Background (CMB) anisotropics 
Bennett et al., 2003; Spergel et al., 2003) should amount to Qy, = (4.6 ± 0.2) 
% of the total cosmological mass/energy density fl^, + ^dm + ^de- While this 
number is in good agreement with the number of baryons 'counted' in the H I 
Lya Forest at z > 2, Qb > 3.5 % (Ranch, 1998; Weinberg et al., 1997), a 
large discrepancy is found in the Local Universe, at z < 2. The total number 
of baryons seen at low redshift in virialized matter (i.e. stars, neutral H and 
He, molecular H, and X-ray emitting plasma in clusters of galaxies: Fukugita, 
2003) and not yet viriahzed matter (residual H I Lya Forest - Penton, Stocke 
& ShuU, 2004 -, and both photoionized and coUisionally ionized OVI in the 
Intergalactic Medium - Tripp, Savage & Jenkins, 2000; Savage et al., 2002) 
amount to only ilb{z < 2) = (2.5 ± 0.3) %. About 46 % of the total predicted 
baryons are then missing at 2; < 2; r2b(missing) = (2.11q;4) % (Nicastro et al., 
2004a - N04a -, and references therein). 

Hydrodynamical simulations for the formation of large-scale structures in the 
Universe predict that these baryons hide in a filamentary web of tenuous 
(nb 10"^ - 5 X 10"^ cm"^, i.e. S = n^,/ < nb >^ 5 - 100 ^ ) mat- 
ter at warm- hot temperatures (T ~ 10^ — 10^ K) that permeates intergalac- 
tic space and connects already virialized structures (e.g. Hellsten, Gnedin & 
Miralda-Escude, 1998; Cen & Ostriker, 1999; Dave et al., 2002; Fang, Bryan 
& Canizares, 2002). This primordial matter would provide the fuel for the as- 
sembly of galaxies and galaxy groups and clusters, and would be shock-heated 
during the continuous infall process of structure formation. The internal shocks 
ionizes the gas to such a degree that it becomes 'invisible' in infrared, optical 
or UV light, but should shine and absorb Far-UV and X-ray photons. However, 
given the very low-density, extremely high sensitivity, large-field of view, FUV 
or X-ray detectors are needed to image the emission from these filaments (e.g. 
Yoshikawa et al., 2003), and these are not yet available. The only avenue pur- 
suable today, is to detect the WHIM in absorption against bright background 
sources. 

Here we present two LETG spectra of the blazars Mrk 421 and lES 1028-1-511, 
and the first detections of three O VH WHIM filaments. A fourth filament is 
only detected in CVis_^2p in the lES 1028-1-511 spectrum. Sections 4 and 5 
are devoted to the estimates of the number density and the cosmological mass 
density in the OVII WHIM, while in section 6 we discuss future prospects 
for WHIM studies. Throughout the paper we use LETG spectra grouped at a 
resolution of AA = 12.5 mA for fitting purposes, i.e. 4 times better than the 
intrinsic LETG resolution of 50 mA. The plots, instead, show spectra binned 
at 25 mA for Mrk 421 and 12.5 mA for lES 1028+511. Equivalent Width 
sensitivity limits are converted into ion column density limits by using Doppler 
parameters corresponding to thermal broadening of the lines in plasma at 



1 < nb >= 2 X 10-'^(1 + zfin^h^ /0.02) is the average density in the Universe. 
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Fig. 1. 21-22.5 A portion of the Chandra-LETG spectrum of Mrk 421. The dashed 
Une shows the best fitting continuum, while the sohd fine is our best-fitting contin- 
uum plus line model. 

logT(K) = 6.1, the central temperature of the WHIM distribution (e.g. Dave et 
al, 2002). We adopt Hq = 72 km s"^ Mpc"^ (Freedman et al, 2001, Bennett et 
al., 2003). Errors are quoted at la for 1 interesting parameter, unless otherwise 
stated. 



2 LETG Spectra of Mrk 421 and lES 1028+511 



Mrk 421 is a blazar at z = 0.03. We observed Mrk 421 with the Chandra 
LETG during two very luminous outbursts, on 2002, October, 26-27 and 2003, 
July 1-2, under our TOO programs to observe blazars in outburst phases. 
These two observations lasted ~ 100 ks each and caught the source at historical 
maxima, 60 and 40 mCrab in the 0.5-2 keV band, allowing us to collect a total 
of ~ 5300 Counts Per Resolution Element (CPREs) at 21 A (Nicastro et al., 
2004b: N04b). This S/N was enough to detect OVII columns of Nqvu^SxIO^^ 
cm~^ at > 3(7. We detected 24 absorption lines in the LETG spectrum of 
Mrk 421, between 10 and 50 A (see No4b for a full, detailed spectral analysis 
of both continuum and lines). While the majority of these lines are imprinted 
by both neutral and highly ionized systems at velocity consistent with zero 
(ISM gas, and either Local Group WHIM - Nicastro et al, 2002; Nicastro et 
al., 2003; Williams et al., 2004 - or an extended Galactic Corona - Sembach et 
al., 2003; Savage et al., 2004), 9 lines are instead identified with two intervening 
WHIM systems at z = (O.OlliO.OOl) and z = (0.027±0.001). Figure 1 shows 
the 21-22.5 A portion of the LETG spectrum of Mrk 421, in which the two 
OVIIi^5_^2p lines from the two WHIM systems, together with the z = lines 
of OVI and OVII (and possibly a fourth O VIIi<j^2p transition at z = 0.033, 
where also a HI Lyct is seen - Dull et al., 1996), are present. 
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Fig. 2. Two portions of the HRC-LETG spectrum of lES 1028+511: (a) 28-30 A 
and (b) 48-50 A. The solid hnes show our best fit models (the continuum is fitted 
over the broad range 10-60 A) . Two hnes from O VII (top panel) and C V (bottom 
panel), from two intervening WHIM systems, are detected in these two portions of 
the lES 1028+511 spectrum. 

lES 1028+511 is a blazar ai z = 0.361. It was observed under our Chandra 
TOO program, with the Chandra LETG, following a short outburst phase 
(3-days), on 2004, March 11. Although lES 1028+511 was back close to its 
'quiescent' level when Chandra pointed at it, lES 1028+511 turned out to be a 
bright soft X-ray source for such a high redshift (-Fo.5_2keV = 1.1 mCrab). The 
150 ks HRC-LETG spectrum of lES 1028+511 contains 80-100 CPREs in the 
21 < A(A) < 50 interval, enough to detect at 3a Nqvh > 1.4 x 10^^(1 + z)'"^ 
cm~^ at 2; < 0.361, and 4.1 times lower CV columns, Ncv > 3.4x 10^^(1 + ;^)"^ 
cm~^ at 2; < 0.269 (at A > 51 A i.e. z > 0.269 for CVi5^2p, the spectrum is 
heavily affected by the presence of the HRC plate gap) . Although the analysis 
of these data is still in progress, we show here some preliminary results. We 
detect two strong lines (3.5(7 and 4.4o-) at A = 28.740 ± 0.038 A and A = 
48.816 ± 0.040 A that we identify as 0VIIi,^2p from a WHIM system at 
z = 0.330 ± 0.002, and C Vi,^2p from a WHIM system at z ^ 0.212 ± 0.001 
(Figure 2 (a) and (b) respectively; Nicastro et al., 2005, in preparation: N05). 
Errors on A and z include a 20 mA (la) calibration uncertainty in the HRC- 
LETG dispersion relationship (see e.g. N04b and references therein). We also 
detect the N VII Lya at the redshift of the O VII system {z = 0.330), but only 
at a significance of 2.6cr. 



3 WHIM Solutions for the O VII Systems 



For the two O VII-N VII absorbers along the line of sight to Mrk 421 (and also 
based on their measured H I and O VI upper hmits, from HST-STIS and FUSE 
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Table 1 

WHIM Solutions for the Three O VII Systems 



Redshift 


T (10^ K) 


Nb « 




[0/H] 


(0.011 ±0.001) 


(0.6-2.5) 


(2.0 ±0.3)^= 


(0.7 ±0.1)^ 


> -1.47 


(0.027 ±0.001) 


(1.1-1.7) 


(2.8 ±0.2)^ 


(0.9 ±0.1)^ 


> -1.32 


(0.330 ± 0.002) 


(1.6-2.5) 


(24 ± 8)"=* 


(7.8 ±2.5)'=' 


N/A 



° In units of io^Qio-[0/H]-i ^^-2_ b ^j^-^g l{]-\o/H]-i^^^^_^^yi (n^.g) is 

the baryon density in units of 10~ cm" "3. " At logT = 6.1. At logT = 6.3. 

respectively - N04b) we find self-consistent ionization balance and metallicity 
solution, that allow us to derive equivalent H column density and tempera- 
ture range (N04a, N04b). These solutions are found by comparing the data 
with hybrid models of tenuous (i.e. ne^lO"^ cm~^) coUisionally ionized gas, 
undergoing the additional contribution of photoionization by the diffuse ex- 
tragalactic UV-X-ray background (plus the proximity effect of the beamed 
radiation from the blazar, in the case of the line of sight to Mrk 421: see N04b 
for details). From these solutions, and assuming homogeneity in the systems 
(which may not be appropriate, given the evidence for the co-existence of 
multi-phase structures for at least one of the two filaments: N04b), we derive 
the depth of these filaments along our line of sight, in terms of the baryon 
volume density in units of nb = 10~^ cm~^. The first two rows of table 1 
summarize the results. 

For the two O VII/CV systems along the line of sight to lES 1028±511, our 
analysis is still in progress. Table 1 (3rd row) contains only rough estimates of 
the equivalent H column density, temperature and depth of the O VII filament, 
derived based on the Nqvii, Nnvii measurements and the Novm 3(7 upper 
hmit. MetaUicity estimates for this object are currently not available. 



4 Number Density of O VII WHIM Filaments 



Figure 4 shows the predicted cumulative number of O VII WHIM system per 
unit redshift as a function of the minimum O VII column density (solid curve; 
Fang, Bryan & Canizarcs, 2002). The lowest O VII column density we detect 
along the line of sight to Mrk 421, is Nq yn = 7 x 10^"^ cm~^, which is about the 
3(j sensitivity of the spectrum of Mrk 421 in the entire z < 0.03 redshift range. 
The redshift of lES 1028+511, instead, is much higher, so we cannot neglect 
the dependence of the Nqvii detection sensitivity threshold on the wavelength 
(see §6). The LETG spectrum of lES 1028±511 has a > 3cr-sensitivity to 
Nqvii = 8.7 x 10^^ cm'^ (the column of the detected O VII line at A = 28.74) 
for 0.269 < z < 0.361, a redshift interval of = 0.092. M z < 0.269 the 
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Fig. 3. Predicted cumulative number of OVII WHIM systems per unit redshift 
as a function of the minimum OVII column density (solid curve; Fang, Bryan &; 
Canizares, 2002). The two points are the cumulative numbers of > 3cr detected 
O VII WHIM filaments per unit redshifts, down to Nq yff = 7 x 10^^ cm~^ (line of 
sight to Mrk 421) and Ng'^ff = 8.7 x 10^^ cm'^ (line of sight to lES 1028+511). 

effective redshift range within which wc can detect Nqvii > 8.7 x 10^^ cm^^ at 
a significance > Sa is Azcs = 0.075 (see N05 for details). So, the total effective 
redshift range is A2;tot = A^^eff + Az = 0.167. 

The cumulative numbers of > 3a detected O VII WHIM filaments per unit 
redshifts, down to N^'^^ff = 7 x 10^^ cm^^ and N^^^ff = 8.7 x 10^^ cm'^, are 
therefore: dAf/dz{> 7 x 10^^) = 70^1^ and dM/dz{> 8.7 x 10^^) = Qtf (Fig. 
3; the large errors are due to the small number statistics: Gehrels, 1986). At 
1.4(7 the two points are inconsistent with a constant value. The best fit slope 
of the power law connecting the two points is a = 1.0to;7. The solid line in 
Figure 3 shows that the model predictions of [14] are consistent (within the 
large errors) in both normalization and slope with the observations. 



5 Cosmological Mass Density of O VII WHIM Filaments 

Based on the equivalent H column density measurements from three OVII 
WHIM filaments (Table 1), we can now estimate the baryonic cosmological 
mass densities l]^^^'^ down to our two values of NQyjf. We obtain r2^^^^(> 
7 X 10^^) = {2.7tll) X 10-[o/^^l-i % (N04a) and Q™^{> 8.7 x lO^^) ^ 
(2.0lt8) xlO"'°/^i-^ % (N05), for the lines of sight to Mrk 421 and lES 1028+511 
respectively. Errors on these estimates include the (dominant) asymmetric 
Poisson errors associated with a small number of events (Gehrels, 1986). One 
can try to reduce these errors combining the two fi^^^'^^ measurements. How- 
ever, to do this, we have to extrapolate the lES 1028+511 measurement, down 
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to the sensitivity of the Mrk 421 spectrum, i.e. down to Nqvii = 7x 10^'^ cm~^. 
First of all we have to derive an equivalent redshift range Az^q, within which 
we would expect to detect a single O VII WHIM filament with Nq yn = 7x 10^^ 
cm-2 along the line of sight to lES 1028+511. We obtain Az^ = 0.023 (see 
N05 for details). Prom Azeq we can derive an equivalent dl^^ distance to 
lES 1028+511 that we will use in the formula for fib- Next we need to eval- 
uate the probability P that this single O VII filament has Nqvii > 8.7 x 10^^ 
cm~^, and use it in the formula for Q]^ to weight the N^^^ measurement (Table 
1). We find P ~ 0.14 (see N05 for details). This gives a combined fib estimate 
qWhim(> 7 X 10^4) = (2.41}:?) X 10-[°/^]-i %, consistent with both model 
predictions and the actual number of missing baryons. 



6 WHIM Detectability 



The WHIM imprints high-ionization metal absorption lines, both in the FUV 
(mainly OVI) and in the X-rays (transitions from He-like and H-like C, N, 
O and Ne, and inner shell transitions from lower ionization ions). While high 
signal-to-noise FUSE FUV spectra of tens of low-redshift AGNs are readily 
available, the observed number density of relatively cool (T 10^ — 10^'^ K) 
shock-heated OVI filaments with Nqvi > 4x 10^^ cm^^ (the average threshold 
sensitivity of high signal to noise FUSE spectra of AGNs) is low: dj\f/ dz ~ 
14lg (Savage et al., 2002). The predicted number density of the hotter (T ~ 
10^'^ — 10^'^) O Vll-dominated filaments, down to the same columns (Nqvii > 
4 X 10^^ cm~^), is a factor of 5 larger, and decreases to that observed for OVI 
at column densities ~ 2 orders of magnitude larger (Nqvii > 4 x lO"*^^ cm~^: 
Fig. 3). The total equivalent H column density (and so the baryon mass) 
imphed by 'pure' OVII filaments (i.e. systems in which the OVI columns 
are below the current sensitivity threshold) is then more than one order of 
magnitude larger than that of 'pure' OVI filaments: the bulk of the WHIM 
can only be detected in the X-ray band. Detecting He-like absorption lines with 
Nion < 2 X 10^^ at A = 22 A (about the restframe wavelength of O VIIis^2p) at 
the spectral resolution of the current X-ray spectrometers (AA = 50 mA for 
both the Chandra Low Energy Transmission Grating - LETG - and the XMM- 
Newton Refiection Grating Spectrometer - RGS), requires X-ray spectra with 
> 600 CPREs in the continuum (Fig. 4, solid line). This requirement relaxes 
significantly at longer wavelength, since Nj^^''^'^ (for given grating resolution 
and CPREs) decreases with A~^. So, for example, at A = 47 A (z = 1.18 for 
the 0VIIi,_,2p, but only z = 0.16 for the CVi,^2p) only - 120 CPREs are 
required to reach the same Scr sensitivity of Nion < 2 x 10^^ (Fig. 4, dashed 
line). Unfortunately, however, the highest quality LETG or RGS spectra of the 
brightest (and very nearby) AGNs contains at most 50 CPREs. Two important 
exceptions are the cases presented here: the LETG spectra of the two blazars 
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LoglN^iL ) 
Fig. 4. CPREs as a function of 3c7 N^^^es^^ ^ 

we use the oscillator strength of the 
0VIIis^2p! / = 0.696). The soUd and dashed lines are the approximate analytical 
relation N^^^es^^ ~ i.i x 10i^iV^(AA(mA))(CPREs)-V2A(A)-^ at 22 A and 47 A, 
respectively. In this formula columns are in cm"^, No- is the number of a (3 in the 
plotted curves), AA is the grating resolution in mA (50 mA for both curves). Data 
points are the measured 3cr ^J^thrmHe-Hke upper limits at 22 A (filled circles) 
and at 47 A (empty circles) from the Chandra spectra of Mrk 421, PKS 2155-304 
and lES 1028+511 and the XMM-Newton spectrum of IH 0414+009. These em- 
pirical data nicely confirm the analytical curves for CFRE^ilOO. At lower CPREs 
data deviate from the ideal analytical curves, due to deviations of the actual sig- 
nal-to-noise ratio in the data from the ideal Poissonian V CPREs. We note that 
about 80 CPREs are needed for a 3o" detection of Nne-Uke = 1-5 X 10^^ cm~^ at 22 
A (lES 1028+511, filled circle), but about 6 times lower columns are 3(T-detectable 
with the same number of CPRE, at 47 A (lES 1028+511, empty circle). 

Mrk 421 {z = 0.03) and lES 1028+511 {z = 0.361), observed by Chandra 
under our TOO programs, while undergoing outbursts. 



7 Future Prospects 



This work only begins the study of the Warm-Hot IGM outside the Local 
group (Nicastro et al., 2002, 2003; Williams et al, 2004), now that the WHIM 
has been detected and theoretical predictions have been verified. 

The large number of ion species showing up in the X-ray spectra imply a rich 
field of investigation. Higher significance detections, along several other lines 
of sight, and at higher redshift, of He- like and H-like transitions from C, N, 
O and Ne, will allow us to measure the relative metallicity ratios at different 
redshifts. This will fundamentally contribute to the assessment of the still 
poor measurements of metal production with cosmic age, and will help refining 
mechanisms for galaxy/ AG N-IGM feedback (e.g. galaxy superwinds, quasar 



8 



winds). Moreover, the potential for the simultaneous detection, in the same 
band, of inner-shell transition from lower ionization species of the same ions, 
will allow us to distinguish between different ionization scenarios, and to assess 
the importance and frequency of multiphase IGM. Spectral resolutions of i? > 
3000 in the soft X-ray band will allow us to separate different kinematics 
components. This has already been proven very useful in the FUV band, 
for example to distinguish between a Local Group WHIM origin and a halo 
origin for the OVI high- velocity absorbers at 2; ~ (Nicastro et al., 2003). 
Reconstructing the distribution functions of several ions in the WHIM (i.e. 
dM I dzdN) , will also probe very useful in distinguishing different cosmologies, 
and possibly tracing the distribution of dark-matter potential wells in the local 
Universe. 

Observing sources with Chandra while in bright outburst phases has proven 
fruitful, allowing the detection of the WHIM for the first time. However, 
Mrk 421 is unique: no other extragalactic sources in the sky (except Gamma- 
Ray Burst X-ray afterglows, e.g. Fiore et al., 2000) reach the level of Mrk 421 
during the observations presented here. While relatively short Chandra and or 
XMM- Newton observations of blazars in outburst will continue to provide a 
'cheap' way to probe new sightlines for WHIM studies, substantial progresses 
in our understanding of the WHIM phenomenon can only be made by exploit- 
ing different observational strategies. 

As shown in Fig. 3, the expected number of OVII WHIM filaments down 
to a given Nqvh threshold, increases hnearly with redshift. About 2 OVII 

filaments with Nqvh > 10^^ are expected up to ^ = 0.5 along a random 
line of sight. Such OVII columns can easily be detected with 100 ks Chandra 
and/or Newton-XMM observations of background sources flaring at a level 
even 10 times lower than that reached by Mrk 421 during the observations 
presented here. 

An important metric for future instrumentation to explore the WHIM, as 
shown in flgure 4, is that the sensitivity of a grating spectrum (i.e. AA = const) 
with constant signal to noise per resolution element, to a given column density, 
increases rapidly with redshift (oc {1 + z)"^). Smaller columns are detectable 
at higher redshift, as beautifully proven by the faint C Vis^2p detection in the 
moderate signal-to-noise Chandra spectrum of lES 1028+511 (about 60 times 
fainter flux than the Chandra Mrk 421 observations presented here: Fig. 2b). 
These are key elements for planning future observational strategies and mission 
designs for WHIM studies. Long integrations, ^0.5 — 1 Ms, with Chandra 
and/or XMM- Newton of high- 2; sources in their quiescent states are low-risk 
and possibly the only way to dramatically increase the number of detections 
and so narrow down the large statistical errors on the fi^^^^'^ estimate. Tripling 
the currently detected number of WHIM systems, by exploiting the above 
observational strategy, is within the reach of current X-ray instruments, and 
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will allow us to reduce the uncertainties on fi^^^'^ down to ~l4o %, the level 
to which the other baryonic components in the local Universe are now known 
(e.g. N04a and reference therein). 

Further advances in this direction will have to wait for the large collecting area 
of the spectrometers on the two planned X-ray missions, Constellation-X ^ 
and XEUS ^ or dedicated small missions, e.g. Pharos (Elvis et al., 2002). With 
~100 times larger effective area and R >1000, Constellation-X will allow us to 
detect WHIM systems as weak as those presented in this work, for dozens of 
different hues of sight, and thicker filaments will be detected towards hundreds 
of different sightlines. 

This will be much more effectively done with gratings than calorimeters, par- 
ticularly if the currently undergoing tests on the off-plane reflection gratings 
configuration for Constellation-X will confirm the recently reached first-order 
peak absolute efficiency of 43 % (Heilmann et al., 2004) ^ . In calorimeters the 
resolution is constant in energy, and so the resolution element in A A degrades 
as oc A^. This dependence exactly cancels out the increase in Nion sensitivity 
given by the formula N'[^^^^{AX)X~^, at AA = const (Fig. 4). Moreover, the 
degradation of calorimeter resolution with wavelength would greatly hamper 
the possibility of resolving WHIM lines (expected at A > 20 A) or detect- 
ing multiphase structures in the IGM (and so prevent tests of the current 
paradigms for the dominant heating mechanism: internal shocks) . This would 
in turn hamper our capabilities of tightly constrain the physical properties of 
the baryons in this diffuse web of WHIM. 

A resolving power R;;5000 is needed to resolve thermally broadened O lines in 
gas at T = 10^ K and so to allow for separation of different components, and 
for a clear physical diagnostics (Elvis et al., 2002), which may be achievable 
with Constellation-X or Pharos. 

Finally we note that it is crucial that the intrinsic noise of the detector con- 
taining the dispersed spectrum be extremely low. Figure 4 shows that the 
lowest CPRE point, from the RGS observation of IH 0414+009, is actually 
closer to the ideal analytical sensitivity curve, CPREs vs N^lJ^^^^j^^,, than the 
one from the HRC-LETG observation of lES 1028+511. The HRC-LETG 
spectrum of lES 1028+511 has almost twice as many CPRE as the RGS 
spectrum of IH 0414+009, but a 5 times larger intrinsic instrumental noise 
(5.2 X 10~^ cts s~^ per resolution element), which greatly penalizes the detec- 
tion of faint columns. High resolving-power/efficiency X-ray gratings dispersed 
onto low-noise detectors are therefore, by far, the most efficient instruments 
to systematically study the WHIM in the near future. 

^ http :/ /constellation . gsfc . nasa. gov /science/index. html 

^ http://www.rssd. esa.int/index.php?project=XEUS 

^ snl.mit.edu/papers/papers/2003/Heilmann/RKH-SPIE-5168.pdf 
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